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Abstract 

We develop the canonical formalism for a system of Af bodies in lineal gravity and 
obtain exact solutions to the equations of motion for N = 2. The determining equation 
of the Hamiltonian is derived in the form of a transcendental equation, which leads to 
the exact Hamiltonian to infinite order of the gravitational coupling constant. In the 
equal mass case explicit expressions of the trajectories of the particles are given as the 
functions of the proper time, which show characteristic features of the motion depend- 
ing on the strength of gravity (mass) and the magnitude and sign of the cosmological 
constant. As expected, we find that a positive cosmological constant has a repulsive 
effect on the motion, while a negative one has an attractive effect. However, some 
surprising features emerge that are absent for vanishing cosmological constant. For a 
certain range of the negative cosmological constant the motion shows a double maxi- 
mum behavior as a combined result of an induced momentum-dependent cosmological 
potential and the gravitational attraction between the particles. For a positive cosmo- 
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logical constant, not only bounded motions but also unbounded ones are realized. The 
change of the metric along the movement of the particles is also exactly derived. 



1 INTRODUCTION 

Solving the motion of a system of N particles interacting through their mutual gravitational 
forces has long been one of the important themes in physics jl] . Though an exact solution 
is known in the N = 2 case in Newtonian theory, in the context of the general theory of 
relativity the motion of the iV bodies cannot be solved exactly due to dissipation of energy in 
the form of gravitational radiation, even when N = 2. Hence analysis of a two body system 
in general relativity (e.g. binary pulsars) necessarily involves resorting to approximation 
methods such as a post-Newtonian expansion 0, ||. 

However in the past decade lower dimensional versions of general relativity ( both in (1+1) 
and (2+1) dimensions) have been extensively investigated from both classical and quantum 
perspectives. Here the reduced dimensionality yields an absence of gravitational radiation. 
Although this desirable physical feature is lost in such models (at least in the vacuum), 
most (if not all) of the remaining conceptual features of relativistic gravity are retained. 
Hence their mathematical simplificity offers the hope of obtaining a deep understanding of 
the nature of gravitation in a wide variety of physical situations. It is with this motivation 
that we consider the iV-body problem in lower dimensional gravity. 

Specifically, we consider the gravitational iV-body problem in two spacetime dimensions. 
Such lineal theories of gravity have found widespread use in other problems in physics. The 
simplest such theory (sometimes referred to as Jackiw-Teitelboim (JT) theory |jj) sets the 
Ricci scalar equal to a constant, with other matter fields evolving in this constant-curvature 
two-dimensional spacetime. Another such theory (sometimes referred to as R = T theory) 
sets the Ricci scalar equal to the trace of the stress-energy of the prescribed matter fields 
and sources - in this manner, matter governs the evolution of spacetime curvature which 
reciprocally governs the evolution of matter || • This theory has a consistent Newtonian limit 
|| (a problematic issue for a generic (1 + l)-dimensional gravity theory 0), and reduces to 
JT theory if the stress-energy is that of a cosmo logical constant. 

The iV-body problem, then, can be formulated in relativistic gravity by taking the matter 
action to be that of N point-particles minimally coupled to gravity. In previous work we 
developed the canonical formalism for this action in R = T lineal gravity J!]] and derived the 
exact Hamiltonian for N = 2 as a solution to a transcendental equation which is valid to 
infinite order in the gravitational coupling constant |J. In the slow motion, weak field limit 
this Hamiltonian coincides with that of Newtonian gravity in (1 + 1) dimensions, and in the 
limit where all bodies are massless, spacetime is flat. 

More recently we have extended this case to include a cosmological constant A, so that 
in the limit where all bodies are massless, spacetime has constant curvature (ie the JT 
theory is obtained), and when A vanishes the situation described in the previous paragraph 
is recovered ||. For A^ = 2, we derived an exact solution for the Hamiltonian as a function of 
the proper separation and the centre-of-inertia momentum of the bodies. In the equal mass 
case an exact solution to the equations of motion for the proper separation of the two point 
masses as a function of their mutual proper time was also obtained. The trajectories showed 
characteristic structures depending on the values of a cosmological constant A. The purpose 
of this paper is to more fully describe these results and to expand upon them. Specifically, 
we generalize our previous formalism with A = J?| to a system of N particles in (1+1) 
dimensional gravity with cosmological constant. When N = 2 we obtain exact solutions for 



the motion of two bodies of unequal (and equal) mass. 

Since the Einstein action is a topological invariant in (1+1) dimensions, we must incor- 
porate a scalar (dilaton) field into the action [|Kj . By a canonical reduction of the action, the 
Hamiltonian is defined as a spatial integral of the second derivative of the dilaton field, which 
is a function of the canonical variables of the particles (coordinates and momenta) and is 
determined from the constraint equations. For a system of two particles an equation which 
determines the Hamiltonian in terms of the remaining degrees of freedom of the system 
is derived from the matching conditions of the solution to the constraint equations. We 
refer to this equation as the determining equation; it is a transcendental equation which 
allows one to determine the Hamiltonian in terms of the centre of inertia momentum and 
proper separation of the bodies. The canonical equations of motion are derived from the 
Hamiltonian. For the equal mass case they can be transformed so that the separation and 
momentum are given by differential equations in terms of the proper time. In this latter 
form they can be solved exactly in terms of hyperbolic and/or trigonometric functions. 

Corresponding to the values of the magnitudes (and signs) of the energy and other param- 
eters (e.g. gravitational coupling constant, masses, cosmological constant) several different 
types of motion are expected in the 2 body case. Broadly speaking, the two particles could 
remain either bounded or unbounded, or else balanced between these two conditions. We 
shall analyze these various states of motion, and discuss the transitions which occur between 
them. We find several surprising situations, including the onset of mutual repulsion for a 
range of values of negative A and the masses, and the diverging separation of the two bodies 
at finite proper time for a range of values of positive A. 

We shall also consider the unequal mass case. In this situation the proper time is no longer 
the same for the two particles, and so a description of the motion requires a more careful 
analysis. We find that we are able to obtain phase space trajectories from the determining 
equation. We also can obtain explicit solutions for the proper separation in terms of a 
transformed time coordinate which reduces to the mutual proper time in the case of equal 
mass. 

In Sec. II we describe the canonical reduction of the theory and define the Hamiltonian 
for the A^-body system. In Sec. Ill we solve the constraint equations for the two-body case 
and get the determining equation of the Hamiltonian, from which the canonical equations 
of motion are explicitly derived. We investigate the motion for A = in Sec. IV, by using 
the exact solutions to the canonical equations. The motion of equal masses for A ^ are 
analyzed in Sec.V where a general discussion on the structure of the determining equation, 
the plots of phase space trajectories, the analysis of the explicit solutions in terms of the 
proper time are developed. We treat the unequal mass case in Sec. VI. Sec. VII contains 
concluding remarks and directions for further work. The solution of the metric tensor, a 
test particle approximation in the small mass limit of one of the particles and the causal 
relationships between particles in unbounded motion are given in Appendices. 



2 CANONIC ALLY REDUCED H AMILTONIAN of N 
PARTICLES 

The action integral for the gravitational field coupled to iV point masses is 
1 = /^[^v^^I^^ + ^V^V^ + ^a} 

+ E / dr a l-m a (-9t,u(x)^^j \ 5 2 (x - z a (r a )) 



1 a "" a 



where \l/ is the dilaton field, A is the cosmological constant, g^ v and g are the metric and its 
determinant, R is the Ricci scalar, and r a is the proper time of a-th particle, respectively, 
with k = 8ttG/c 4 . The symbol V M denotes the covariant derivative associated with g^„. 
The field equations derived from the action ([!]) are 

R-g» u V,V^ = 0, (2) 

-V M tf V„tf - -^„V A W A # + ^,V A V A ^ - VpVvV = kT^ + - 9lu ,A , (3) 

4- {^^f\ - \guU**)^r^r = ° ' (4) 

dr a { dr a J 2 dr a dr a 

where the stress-energy due to the point masses is 

T^ = y2m a dT a —=g^g up —2-—2-5 2 (x - z a {r a )) . (5) 

a J V~9 dr a dr a 

Eq.flU) guarantees the conservation of T^ v . Inserting the trace of Eq.fl3|) into Eq.(||) yields 

R-A = kT^. (6) 

Eqs. (HD and (^|) form a closed sytem of equations for the matter-gravity system. The 
evolution of the dilaton then follows from inserting the solutions to these equations into(Q), 
and then solving for its motion, the traceless part of (|3]) being identities once the other 
equations are satisfied. Alternatively, one can solve the independent parts of equations (@), 
(|3]), and (|j) for the metric, dilaton and matter degrees of freedom, which is the approach 
we shall take in this paper. If the masses of all particles are taken to be zero then the field 
equations reduce to those of constant curvature lineal gravity, or JT theory j|, Q . 

Consider next the transformation of the action (jl]) to canonical form. We decompose the 
scalar curvature in terms of the extrinsic curvature K via 

V^gR = -2d (^K) + 2d 1 (^N l K - r l diN ) , (7) 

where K = (2N 'y)~ 1 (2d 1 N 1 — 7~ 1 A 1 <9i7 — cVy), an d the metric is 

ds 2 = -N$dt 2 + 7 j dx + —dt J , (8) 



so that 7 = gi±,N = (— g 00 )^ 1 / 2 and N\ = gio, and then rewrite the particle action in 
first-order form. After some manipulation we find that the action (|I|) may be rewritten in 
the form 

1 = Jdx 2 lj2P^5(x - z a (t)) + vr 7 + n^ + N R° + N.R 1 ] , (9) 

where ir and IT are conjugate momenta to 7 and \1/, respectively, and 

a V I 

(10) 
R 1 = ^7r--TM' + 27r' + Y^5(x-z a {t)) , (11) 

7 7 a 1 

with the symbols ( ' ) and ( ' ) denoting d and di, respectively. 
The action (H) leads to the system of equations 

* + N — v^tt 2 -— vrn+— -— vl/' 2 - _ 

I 2 v/7 8/c/rY 4^7 k 

a 27 2^ +m 2 

+ N, {—!& + - + £ ~ 2 S(x - z a (t))\ + K-±—V + N[- = , (12) 
I 7 7 a T J 2*^77 7 

7 - iV (2/^77*-- 2k v^I) + iVi--2iV{=0, (13) 

7 

R° = , (14) 

i? 1 = , (15) 

ri + ^(-i^n + — ^-W + -^) = , (16) 

4? + N {2KJyir) ~ N^-V) = , (17) 
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P« + ir : t- + m «- 



fe«f 7 a 2 M + m 2 7 2 ^a <^a 7 



+A ^Sr - (18) 

z a - iV /2 7 + — = . (19) 

EL + m 2 7 



In the equations ( |i"8|) and (|I9|), all metric components (No, N%, 7) are evaluated at the point 
x = z n and 

df _ df(x) 



-a 



dz a dx 

This system of equations can be shown to be equivalent to the set of equations (||), (|3]) and 

Since iVo and Ni are Lagrange multipliers, equations ([TJ]) and (|1"5| ) are constraints; specif- 
ically they are the energy and momentum constraints of the (1 + 1) dimensional gravitational 
system we consider. We may solve them for (^f'/y/^y)' and it' in terms of the dynamical and 
gauge (i.e. co-ordinate) degrees of freedom, since these are the only linear terms in these 
constraints. We identify these coordinate degrees of freedom by writing the generator arising 
from the variation of the action at the boundaries in terms of (^f'/y/^y)' and it' , and then 
finding which quantities serve to fix the frame of the physical space-time coordinates in a 
manner similar to the (3 + l)-dimensional case. 

Carrying out the same procedure as in the A = case || we find that we can consistently 
choose the coordinate conditions 

7=1 and IT = . (20) 

Eliminating the constraints, the action reduces to 

1 = j dx 2 r£ Va z a 5(x- Za )-n\ , (21) 

and the reduced Hamiltonian for the iV-body system is 

H = [ dxTi = -- [ dxAV , (22) 



H 

where \1/ is a function of z a and p a , determined by solving the constraints which are under 
the coordinate conditions ( p0|) 



A*- -(V) 2 + ^ ~2+ K £ M + ml5(x-z a ) = 0, (23) 

27i' + J2PaS(x-z a ) = 0. (24) 



The expression for the Hamiltonian fl22|) is analogous to the reduced Hamiltonian in 
(3 + 1) dimensional general relativity. In (1 + 1) dimensions it is determined by the dilaton 
field at spatial infinity. The consistency of this canonical reduction may be demonstrated in 
a manner analogous to that employed in the A = case: namely the canonical equations of 
motion derived from the reduced Hamiltonian (^) are identical with the equations (0) and 

The methodology at this point is then as follows. First we must solve ( |2"3T ) and ( P4Tj for 
\1/ and 7r in terms of the p a and the z a , consistently matching solutions across the boundaries 
of the particles. Then we compute from ( ^21) the Hamiltonian in terms of the independent 
momenta and coordinates of the particles. This expression is sufficient to obtain the phase- 
space trajectories for a given set of initial conditions. Finally we solve equations (p~2| — P~9D to 
obtain a complete solution for the A^ body system. Throughout the remainder of this paper 
we shall consider only N = 2, i.e. 2-body dynamics. 
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3 SOLUTION TO THE CONSTRAINT EQUATIONS 
AND THE HAMILTONIAN FOR A SYSTEM OF 
TWO PARTICLES 

The standard approach for investigating the dynamics of particles is to derive an explicit 
expression of the Hamiltonian, from which the equations of motion and the solution of 
trajectories are obtained. In this section we explain how to derive the Hamiltonian from the 
solution to the constraint equations ( p3|) and ( p4|) and get the explicit Hamiltonian for two 
particles in a spacetime with a cosmological constant. 
We first express the equations fl2"B| ) and ( PH ) as 

1 



A* = - m 2 - k 2 (xT + 2 A - K E pi + ml6(x - z a ) , (25) 

A X =-Wpa^-^a), (26) 

^ a 



where we set x' — it. Rewriting (^5|) as 

(1 + £)A¥ = Ia(M/ 2 - 4kV) + -A + k 2 x A x -kJ2 \jvl + ml5{x - z a ) , (27) 

4 8 I a 

and using (|26|), we can rewrite ( |22"D as 



'? 1 + \9(z a ) 2^1 + \9{Za) 



1 /■ . 1 



dx f— -A(^ 2 -4k 2 x 2 + 2Ax 2 ) , (2? 

1 + |#(x) V ' ' v 



which can also be obtained by inserting (|25| ) and (|26| ) into ( [22] ) and iterating by partial 
integration (assuming convergence). We shall refer to this formula later when we consider 
boundary conditions. 
Defining <fi by 

tf = -41og|0|, (29) 

the constraints ( ^5|) and (|26|) for a two-particle system become 



A0-i{/, 2 (x') 2 -^A}0 = £{Vrf + m!*(siW*-*i) 



+ v /p| + m|0(2; 2 )5(x-2 2 )) , (30) 



A X = -g {pi<J(a;-2i) +p 2 5(x-z 2 )} • (31) 



The general solution to (|3~l"D is 



X = -j {pi | x - zi | +p 2 \x - z 2 \} - tXx + eC x . (32) 



The factor e (e 2 = 1) has been introduced in the constants X and C x so that the T-inversion 
(time reversal) properties of x are explicitly manifest. By definition, e changes sign under 
time reversal and so, therefore, does \. 

Consider first the case z 2 < Z\, for which we may divide space into three regions: Z\ < x 
((+) region), z 2 < x < Z\ ((0) region) and x < z 2 ((-) region). In each region x' is constant: 



-eX-U Pl +p 2 ) 
y \ -eX + U Pl -p 2 ) 

-tX + \( Pl +p 2 ) 



(+) region, 
(0) region, 
(-) region . 



(33) 



General solutions to the homogeneous equation A0 — ~ I k 2 (x') — |A J- = in each region 

(34) 



arc 



io(x) 

Ux) 



■■ A + e 1 2 K + x + B + e~2 K+x , 
A e3 KoX + B e-s Kox , 

. A _ e \K.x + Be -\K.x 



where 



K + 



l^(x + l( Pl +p 2 )) -|A 



rf» (X-U Pl 



V2. 



|A 



tflX-hfa 



V2, 



|A 



(+) region 
(0) region , 
(-) region . 



(35) 



For these solutions to be the actual solutions to Eq.(p0|) with delta function source terms, 
they must satisfy the following matching conditions at x = z±, z 2 : 



<f>+(zi) = M z i) = 000 > 

4>-{z 2 ) = <Pq{z 2 ) =4>{z 2 ) , 

<j)' + (zi) - 0qOi) = -yjpl + m\<\>(z x ) 
fiofa) - <j>'-(z 2 ) = -\j P \ + m 2 2 (j)(z 2 ) 



The conditions fl36a| ) and ( |36c| ) lead to 

e? K+Zl A + + e- 1 * K+Zl B + = e^ K ° Zl A + e^ K ° Zl B 



(36a) 
(36b) 

(36c) 
(36d) 



(37) 



yielding 



.4, 



,JK + zx A 



e 2 



\K + z\-q 



KJ P i + mi + 2K + 2K + L 



KJ P i + mj + 2K i 



2K, 



-e* 



Kqzi 



A 



kJ P { + m{- 2K _i 
H — r^ e 2 



2fG 



-Kqz\ r> 

-DO ) 



4K, 



i(^o-^+)^i^ o 



(38) 



K\p{ + m(- 2K + 2K + i 

_v^i 1 — u i e - 5 (^o+x + ) 215 39a 

47T , 



Hj 



«VPi + m i + 2K o ~ ZK + i 



4K 4 



-e 2 



(K- +i^+)2i^4 



/cy'pf + m( - 2K - 2K^ 
AK~ 



e -±(K -K + )z 1Bo _ (gg^ 



Similarly from ( |36b| ) and (36d) we have 

e^- Z2 A_ + e-^-^B. = e? KoZ2 A + e~^ K ° Z2 B 



(40) 



e |ir_ Z2 



A_ + e~ 



4^-22 



5_ 



ka/P^ + m 2 , - 2K i 



2K. 



et Xo22 A 



+ 



KJpi + mi + 2K _i 



Ko«2 



2K_ 



e -2^Z2 Bi 



and then 
.4_ 



KJP2 + ^2 - 2i^o - 2K_ i 



4K_ 



i(K -^-)^2^ 



nJpj + mi + 2K - 2K_ x 



(41) 



AK. 



e --(Ko+K-)z 2BQ ^ ^ 42f 



5_ 



«a/P^ + m 2 2 - 2K + 2K_ L 



AK. 



±(K +K_)z 2 a 



K\pl + m 2 2 + 2K Q + 2K_ Uv v . 
+ -^ TIT " e-^ K °- K - )z2 Bv . (42b) 



Since the magnitudes of both and \ increase with increasing |x|, it is necessary to impose 
a boundary condition which guarantees that the surface terms which arise in transforming the 
action vanish and simultaneously preserves the fmiteness of the Hamiltonian. A consideration 
of (|2"8|) implies that we may choose the boundary condition 



^ 2 - Ak 2 x 2 + 2Ax 2 = C±x for (+) and (-) regions (43) 

with C± being constants to be determined. This boundary condition means 



A_ = B + = , 

{2K + x + 41og|A+|} 2 - Ak 2 j- 

+ 2Ax 2 = C + x , 
{2K_x - 41og|B_|} 2 - An 2 j- 

+ 2Ax 2 = C_x . 



eX + -{ Pl +p 2 ) 



e x - ^(Pi+P2) 



1 l 2 

x + eC x + -(pizi + p 2 ^) } 

1 v 

x + eC x - -{piZt + P2Z2) > 



(44) 



(45) 



(46) 



It may seem that instead of ( [44"D we could have made the alternate choices (A + = 0, B_ = 
0),(A+ = 0,A_ = 0) or (B + = 0,5_ = 0). However the definitions Q35|) imply that 
K± are positive quantities, which in turn leads to a negative Hamiltonian for the choice 
(A+ = 0, B_ = 0). For the choices (A+ = 0, A_ = 0) and (B+ = 0, 5_ = 0) the Hamiltonian 
identically vanishes. 



The terms quadratic in x from (|45| ) and (051) merely recapitulate the definitions of -KT±. 
Equating terms linear in x yield the relations 



l§K+\og\A + \ + 8k 2 
-16AT_log|5_| +8k 2 



eX + -(pi+p 2 ) 



eCx + T^i+Pa^) 

eC x - 4(^1^1+^222) 



<?4 



Equating the constant terms of ((45]) and (Eg) leads to 



16(log|A, 



4/c 2 



16(log|£_|) 2 -4K 2 



eC< x + 4(^1^1+^222) 
eC x ~ 4(^121+^222) 



i2 








(47) 
(48) 

(49) 
(50) 



We choose the solutions 



log|^ 



log|£- 



e 

4 l 



C x + -(piz 1 +p 2 z 2 ) 



(51) 



C x - -{pizi + p 2 z 2 ) 



Before proceeding, we add a remark to ([5*TD. In solving ( (4D|) and (|5T)|), there are actually four 
combinations (=F, ±) of sign choices for log|^4 + | and log|-B_|. However the choices (+, +) 
and (— , — ) do not lead to any relations among the coefficients and gives us an unphysical 
Hamiltonian. The choices (+, — ) and (— , +) lead to identical physical results once the signs 
of the momenta Pi and the coefficient C x are reversed. 
The condition (HIT) leads to 



and 



A 

Bo 



Bo 



KJpf + m(- 2K - 2K A 



-K q z-l 



nJpi + m( + 2K - 2K 4 



nJpi + m 2 2 + 2K - 2K_ 



-K z 2 



KJpi + mi- 2Ko - 2K. 



(52) 



(53) 



From (E2T) and (153T) we obtain 



K 



p\ + m\- 2K - 2K A 



K 



pl + m 2 2 - 2Ko - 2K. 



KJpf + m{ + 2Ko - 2K + nJp 2 2 + m 2 2 + 2K - IK- e 



,Ko{zx-z%) 



(54) 
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which we shall refer to as the determining equation for X. The \1/ fields in (±) regions are 

\£ + (x) = -4log\A+\ - 2K+x , 

(55) 
*_(ar) = -A\og\B„\+2K-X , 

and the Hamiltonian is 

H = -I/ rfx A* = -V]!°oc 

= 2(AV + /f - ) . (56) 

Once the solution for X is obtained from fl5"J|), the Hamiltonian is explicitly determined from 
(|56"D in terms of the degrees of freedom of the system (i.e. the coordinates and momenta of 
the particles). 

^From (|39a| ), ( |42b|) , (|52"D and (0), A + and S_ are expressed in terms of A as 

A + = iKo , e^ K °- K +^A , 



2K + 2K + - K^p{ + m{ 

(57) 

B _ = 4K e ±(K 0+ K-) Z2A 



2K - 2K_ + K^Jpi + mi 

and from ( |5"TD and ( p7| ) the coefficients A and C x (and hence A + , B_ and -Bo) are also 
determined 



1 (2K + 2K + -K^pl + m 2 1 )(K^pl + ml + 2K -2K + ) 1 

<?x = 7T l °g I l + 7r( K + z l + K - Z V ' 

2k {2K + 2K_-K^jpl + m 2 2 )(K^pl + ml + 2K -2K_) 2k 

(58 



1 2tf + 2iC - Kdpi + m 2 2 1 «; e 1 Ke 

log ^ = -log V -(KoH-#+ + V Pi)*iH-t(#o-#-- V P2>2 

2 «Jpi + ™i + 2if - 2K+ 4 2 4 2 



(59) 



1 2if + 2K + - K\ pi +mj 1 ke 1 ke 
logS_ = -log ° t — --^0-^ + ^1 + 7(^0+^-^2 

2 KJpi + m| + 2K - 2K_ 4 2 4 2 



(60) 

1 (2K + 2K + - K ^/pf + mJ)( K Jrt + ml + 2K - 2iT_) 1 Ke 

loglAol = -log * ^f + l(K+ -K --p 1 )z 1 

-l(K + K_ + ^p 2 )z 2 , (61) 



1 (2K + 2K_-«: v /p| + mi)(«: v /pf + mf + 2Ko-2ir + ) 1 Ke 

log|S | = -log V __| + l(Ko + K+ ._ pi)zi 

-±(K--K + ^- P2 )z 2 . (62) 
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The parameters C± are determined from (|47|) and ( 48 ) . From fl52|) , (|53|) and fl57|) it is evident 
that an overall common sign of A + , 5_, A and i?o has no physical meaning, and so we can 
choose all these coefficients to be positive. 

The previous expressions are somewhat cumbersome. We can express them more com- 
pactly by making use of the following notation: 



K x 



2K Q + 2K_ - nJpi + m\ 



K 2 


= 


2K + 2K + - «VPi + mj 


K 01 


= 


K€ 
Kq- K + + —p! , 


K m 


= 


Kt 

K -K_- -p 2 , 


Mi 


= 


K \ 
K \ 


/pi + mi + 2K 


-2K + 
-2K_ 


M 2 


Ipl + m 2 2 + 2K 


Y + 


= 


K 


x+ e -( Pl +p 2 ) 


? 


Y 


= 


K 


X--{pi-p 2 ) 


? 


Y_ 


= 


K 


x- e -( Pl +p 2 ) 





The coefficients A+, B_, Aq and Bq can then be rewritten as 

K,\ 1/2 



B_ 

A, 
B 



Ms 



e —\(K xzi-Ko2Z2)—\K + zi 



e 4 



\{K 01 z 1 -K 02 z 2 )+\K_z 2 



(K 2 M 2 ) 1/2 -\(K 01 zi-K 2Z 2 )-hK Z2 
4K 



-e 4V 



and the solution for <p is then 

( Kl \ 2 ^\{K al z 1 -K 02 z 2 )+\K+{x-z 1 ) 



\M X 

1 

K 2 



(63) 



(64) 



e 4 



1 { K 01Z1 -K 02 z 2 ) ^ KiMi Y/2 e -\K (x~z 1 ) + ^ K2M ^in e \K (x-z 2 )^ ^ ^ 



M, 



— \ (K (n z 1 -K 02 z 2 )-±K- (x-z 2 ) 



Repeating the analysis for z\ < z 2 yields a similar solution with pi — > —pi. Hence 
the full solution is obtained from the preceding expressions by replacing pi and Z\ — z 2 
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by Pi = Pi sgn^x — z 2 ) and \z\ — z 2 \, respectively. The determining equation (|54]) of the 
Hamiltonian is then expressed as 

K X K 2 = M X M 2 e Ko]Z1 -^ . (66) 



or 



4AT 2 + [nyjpi + ml - 2K + ] [K\fp~l + m\ - 2KJ\\ tanh ( -K Q \zi - z 2 \ 
= -2K ([Ky/rf + ml - 2K+] + [n^/pj+m^ - 2Kj\ , 



(67) 



where the momentum pi is replaced by p^. 

For the expression ( p8|) to have a definite meaning as the Hamiltonian, K± should be 
real. This imposes the restriction H 2 + 8A/K 2 > 16 (pi +P2) 2 - However Kq need not be real. 
If A takes a sufficiently large positive value Kq will be imaginary and the above analysis 
must be repeated. In the (0) region the soluton to the equation (|30|) becomes 



b (x) = A s sm-K x + A c cos-K f 



0% , 



(6f 



where 



K n 



-tKa 



-A-tf\X--(p 1 -p 2 ] 



(69) 



Under the same matching conditions ( |36a[ - |36d| ) and the boundary condition (|43| ) we get, 
instead of fl6~7|), a new determining equation for the Hamiltonian 



i h'C, ■ W\/i>l + ml - 2K+\[H\/i>i -\ nil - 2A'_] ) tan ( -K Q \z x - z 2 \ 



2K [kJpI + ml- 2K + ] + [nJpl + m\ - 2K_ 



(70) 



which is just the equation derived from (|6T| ) by formally replacing K$ with iKq. Similarly, 
the solution for for imaginary Kq is also identical with that obtained from fl65|) by the 
same replacement. Hence equation (^7|) is valid for all values of K , and may be regarded as 
a transcendental equation which determines H as a function of the independent coordinates 
and momenta of the system. 

We have previously shown that in the case of zero cosmological constant the solution for 
H can be expressed in terms of the Lambert W function. In this more general case with 
A /0 the solution for H from ( |B"7|) cannot expressed in terms of known functions. Rather 
we must regard H as being implicitly determined in terms of the coordinates and momenta 
via (|3). 

Finally, the components of the metric may be computed from the equations (|T2D , (13), 
( JIB]) and ( |TTD under the coordinate conditions fl20"|). The derivation and the explicit solutions 
of the metric are given in Appendix A. 
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The canonical equations for the 2-body system can be derived by differentiating the 



determining equation 
z\ this yields 

Pi = 



with respect to the variables Zi and pi . For the variables pi and 



OH 



_2 fdK ± dK. 

dz\ k \ dz\ dz\ 

2 (Y+_ Y_} KpK x K 2 



y± 

K 4 



~K~) ~dz[ 



K 



AT_ 



J 



(71) 



z\ 




y± 



(72) 



where 



J 



K K, 



-2 ^ 



K 






\ 1 



Y 



o 



Y_ 
Uo 



#-> 



Similarly, for particle 2 the equations are 

2 /X Y_ ^ K0K1K2 
K. 



Y_ 



M, 



K X K 2 - ^K x K 2 { Zl 



*2 



P2 
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K 



—e- 



Y_ 



+ 



K -J 

J\K 4 



J 
Y. 

k1 



k k 2 

Mo 



Pi 



Y_ 



(73) 

(74) 
(75) 



Ip2 + ™2 

It is straightforward to show that these canonical equations guarantee the conservation of 
the Hamiltonian (i.e. H = 0) and the total momentum p\ + p 2 (i.e. pi + p 2 = 0) . 

Alternatively, the equations of motion ([18]) and ([19]) derived from the action (H) become 



Pa 



under the coordinate conditions 
given in the Appendix A into 



dN 
dz a 



p 2 a + ml + 



dz. 



-Pa 



No 



Pa 



pl + m 2 a 



-JVi 



(76) 

(77) 



Insertion of the solutions of the metric components 
and ( |77|) reproduces the canonical equations of motion 



(|7l"D , (|72|) , ( f74|) and (]75|) when the partial derivatives at zl,z2 are defined by 



dN 0A 
dzi 



dN 0<1 



dx 



+ 



«9AV 



dx 



(7* 



Thus consistency between the geodesic equations and the canonical equations of motion is 
explicitly verified, while the formal proof of the consistency in the case of A = can be 
easily generalized to the A^O case. 
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EXACT SOLUTIONS OF THE TRAJECTORIES 
IN THE A = CASE 



In a previous paper [§] we showed that in the A = case the determining equation ( |66D 
can be solved explicitly and the Hamiltonian for the equal mass is expressed in the center of 
inertia (C.I.) system p\ = —p 2 = p as 



W 



"f (\ r \Vp 2 + m2 — epr)exp (| (\r\\fp 2 + m l — epr 



H = \ lp 2 + m 2 + ep sgn(r) — 8 — 

v ' K\r\ 

(79) 

where W(x) is the Lambert W function defined via 

y . e y = x => y = W{x) . (80) 



The W(x) has two real branches Wo and W-\ for real x fTT |. 

The Hamiltonian (|79|) is exact to infinite order in the gravitational constant and for 
arbitrary values of m and p. We can view the whole structure of the theory from the weak 
field to the strong field limits. By setting H = Hq we can draw a phase space trajectory in 
(r, p) space. This phase space trajectory should be, as a matter of course, obtainable directly 
from the solution r(t),p(t) to the canonical equations by eliminating the time variable t. 
Indeed, this can be verified by numerically solving the equations 

*(H-2ep)(H-ep-Vp r T^) 
P = ~ 2-f(H-ep-^ + m 2 ) Sgn(r) ' (81) 

[ 2 — yk™ ~ e P ~~ vP + m ) vP + m ) 

which are obtained in the case of A = from (fTT|), (|7^), ([73) and (|75|). However for certain 
values of the parameters, superficial singularities appear in r(t) and p(t) due to the zero 
points of the denominator 1 2 — nr(H — ep — \/p 2 + m 2 )/4>. These singularities correspond 
to W(x) = —1 representing the transit point between two branches W and W-\. In a 
spacetime description the singularities are coordinate singularities and are a consequence of 
t being a coordinate time. 

We can deal with this problem by describing the trajectories of the particles in terms 
of some invariant parameter. The natural candidate is the proper time r a of each particle. 
From the metric components given in Appendix A and the canonical equations (|77|) , the 
proper time is 

dr 2 a = dt 2 {N (z a ) 2 - (JVi(* a ) + z a f) 

pi + m i 

For the equal mass case it is common for both particles 



dt 2 N (z a ) 2 —^- (a = 1,2) . (83) 



(H — 2ep)m , ,„.. 

dr = dn = dr 2 = -. -A: FJ dt , (84) 



{2 - f(H - 2ep)} (Vp 2 + m 2 - ep)V^+ 
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m 2 



and the canonical equations ( JSID and (82) become 
lp 2 + m? < H(Jp 2 + m 2 — ep) — 

ep) 



K 

Am 



dp 

dr 

— = —{\Jp 2 + m 2 



m > sgn{r) 



2e 
m 



K I ■ 

4 



(H — ep — \/p 2 + m 2 ) y/p 2 + m 2 



H -2ep 



(85) 



1 } sgn(r) .(86) 



Remarkably, the equations (p5|) and ( |8"q ) have an exact solution. We can obtain it in the 
following way. First we solve Eq . (JS"5|) for p(r). From this r(r) can be extracted by directly 
solving fl3"E| ) after substituting the solution for p or by solving ( |5"6J ) for r. This yields an 
exact expression for the proper separation r of the two bodies as a function of their mutual 
proper time. Note that since 7 = 1, at a fixed time dt = (and hence a fixed dr = 0), the 
separation r = z± — z 2 is the proper distance between the two particles. 
In the r > region the solution is 



p(r) 
r{r) 



t ( /o(T) " am) 

16 
«{if-m(/o-£)} 



57) 



tanh 



H — m 



(/» + *) 



H-mU 



h 



with 



/»(r) = ^ 

m 



Pq + m 2 - ep 



?9) 



where po is the initial momentum at r = r . In the r < region the solution is 



p(r) 
r(r) 



em 



/o(r) 




-16 



«{ff- m (/ -^)} 



tanh 



-1 



H 



m 



f° + is) 



H 



(f°~~k 



ml Tn — 



with 



Mr) = - 

m 



Po + m- 









e 4 



Po + m z + ep 



(90) 
(91) 

(92) 



In Fig.l and 2 we show the typical plots of r(r) and phase space trajectories for various 
values of H with m = 0.5. (In all plots in this paper we choose k = 1.) In Fig.l it is 
seen that as H increases, not only the amplitude and the period of the bounded motion 
become large, but also the shape of r(r) deforms from the parabolic (Newtonian) shape. 
This deformation is also present in the phase space trajectories shown in Fig. 2. At higher 
energy the trajectories become more and more S shaped. This is due to the fact that the 
trajectory smoothly moves over the two branches W and W-\ [|]. 

We see in Fig. 3 the striking distinction between the separations of the particles in the 
non-relativistic and relativistic cases once the value of the energy becomes large relative to 
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the mass (here H = 25 and m = 0.5). The maximal separation of the particles is much 
smaller than in its Newtonian counterpart (the dashed curve), and is achieved far more 
quickly. After maximal separation, the particles move toward each other at a slower velocity 
until they are very close together. At this point (less than 10% of their maximal separation), 
they accelerate toward the same point, after which the motion repeats with the particles 
interchanged. 



Fig.l 
The exact r vs r curves in the case of A = for m = 0.5 and four different values of H . 



Fig.2 
Phase space trajectories corresponding to the r(r) curves in Fig.l. 



Fig.3 

The exact r vs r curve for H$ = 25 with m = 0.5 and the Newtonian curve for the same H$ 



5 EXACT SOLUTIONS OF THE TRAJECTORIES 
FOR THE EQUAL MASS IN THE A ^ CASE 

5.1 General Discussion 

In this section we consider a system of two particles with equal mass for the A^O case. In 



the C.I. system, depending upon the sign of (\/H 2 + 8A/k 2 — 2ep) 2 — 8A/k 2 the determining 
equations QS7| ) and ( |T0"D become 



2 , d 2n i K 



(Ja + B 2 ) tanh i-J A \r\) = 2J A B , (93) 

and 

(J A 2 - B 2 ) tan (Zj A |r|) = -2J A B , (94) 
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respectively, where 



Jk 



Jk 



\ 



8A 

~~2 



'H 2 + ^- 2ep I - 



K 



8A 



\ 



8A / / „ 8A 



K- 



Equation 



B = H - 2Jp 2 + m? . 
may be further divided into two types: 

B 



tanh ( — JaI)" 



Ja 



(tanh-type A) 



Gi- 



ft \ /7a 

tanh ( —Jx\r\\ = — (tanh-type B) . 



(95) 



(96) 



(97) 



When A = the tanh-type B equation is excluded because J\/B exceeds 1. When a 
cosmological constant is introduced, this equation has solutions in some restricted range of 
the parameters. 

Likewise, eq.(|94]) may also be divided into 



. K ~ . 

tan( j^JaFI 



B 

-=- (tan-type A) 

Jk 



or 



tan ( —Jh\r\\ = -^ (tan-type B) . 



(9S 



(99) 



For all of these four types of determining equations the canonical equations of motion 
are identical: 



P 



kJKJI - B 2 



16C 



sgn(rj 



2eWl 



8A / j£\ zj£ p 

^ I 1 " ^J sgn(r) + -^^^^ • 



(100) 
(101) 



where 
C = 



1 + 



8A 

K 2 H 2 



1 + 



K 



8A 9 / / n 8A \ / k , , 



2^2 ^A 



fi'" 



2ep) (SH--(^-^)r 



(102) 



5.2 Phase-space Trajectories 



For a given value of A the equations (|93|) or (p4j) describe the surface in (r, p, H) space of 
all allowed phase-space trajectories, from which the (r,p) trajectory is obtained by setting 
H = H . 
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Fig.4 shows phase-space plots for a small H = 2.2m and three different values of A under 
identical initial conditions. First we note that the trajectories of the relativistic motion (solid 
curve) are are slightly distorted compared to the Newtonian motion (dashed curve). Second, 
a trajectory with A > (dash-dot curve) is expanded, reflecting the repulsive effect of 
the cosmological constant, while a trajectory of A < (dotted curve) shrinks, due to the 
additional attractive effect. The relativistic plots in Fig.4 correspond to the choice of e = 1. 
Plots for e = — 1 (the time- reversed solutions) are obtained by reflection in the p = axis. 
The phase space plots for a large H Q = 8m are shown in Fig. 5. The trajectories for A > 
become extremely S"-shaped, while for A < the trajectory is still a distorted oval due to 
the attractive effect of A. The effects of the cosmological constant (A > repulsive; A < 
attractive) are precisely analyzed in terms of the exact r(r) plots in the next subsection. 



Fig.4 



Non-relativistic (Newtonian) and relativistic trajectories form = 0.5, Hq = 1.1 and three 
different values of A. The undistorted oval (which is really a pair of parabolas intersecting 

at the r = axis) is the Newtonian trajectory. 



Fig. 5 

Relativistic trajectories for m = 0.5, H = 4 and three different values of A. 

5.3 Explicit Solutions 

The phase-space trajectories discussed in the previous subsection can be obtained from the 
solution to the canonical equations ( |100|) and ( |101| ). For the equal mass case there is a 
common proper time for both particles 

dr = dn = dr 2 = - ™ ,. ^-dt , (103) 

Vp + m z C 



via which the canonical equations ( |100| ) and ( |101| ) may be expressed in the form 



— sgn(r) , (104) 



(105) 
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First we solve Eq. (|104|) for p and then (p6|)- (|99|) for r. In the r > region Eq. (|104|) leads 



to 



p 

PO 



dp 



kH 



{VW+^-^i + ^p-f }Vp f + 



(It 



77T 



provided the condition 



8A 

1 + ^ a0 



Am J T 




kH 

4^ (T - ro) ' 


(106) 




(107) 



is satified. Hence for A < the motion is allowed as long as H satisfies 



H>\ - 



8A 



K 



2 ' 



(10* 



We perfom the integration of the LHS of (|106|) for three separate cases which depend on the 
value of A relative to m and H. The solution p(r) is 



with 



/(' 



*>-TK7k 






H < 



m-ai^(T-T ) 



^+VW) 



1+V^TV 



_cr + -^-\/ — 7m tan 



-Tm(T-To) 



where 



o, — 1 _i_ 8A 



7r> 



1 + 



8A 






n- m 



V 



II 



/111 



• 



with po being the initial momentum at r = To- 
Similarly the solution in r < region is 



*>~T*>-W5 



with 



7 \ g^m. /xt^vZ — ITT 

i+V7^+(v^r-i) v e -3-v^r(— o) 



70 



l + V / 7H 



7m > , 
7m = , 

7m < , 



a = {l + ^)Upl + m^-ep )-^ 



7m > 
7m = 



(109) 



;no) 



:m: 



(112) 



(113) 



r(l+V75") 



1 + V — 7m T- 



"FT^ 



-7m(T-T ) 

-lm(T-ij) 



7m < 
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where 



o- = (1 + ^/ln)Nvl + m 2 + ep ) 



m 
If 



V 



a ~ 7TV7-. 



H V ' m 



(114) 



The solution for r{r) for each of the determining equations (|96|) - (|99|) is obtained as 
follows 
tanh-type A: 



r(r) 



1(5 



( V , K 2 -H' 2 +8A-mK(/(r)- jK)) -8A 



tanh 



^(g-ml/M+^yl) 



-16 



tanh * 



v / K 2 H 2 +8A-mK(/(r)-ji-)j -8A 



(V« 2 -ff 2 +8A-mK(/(T)- jfa)) -8A 



tanh-type B: 




1(5 



yJ(Vii"m-^K-mK{f(r)- ? ^ y )) 2 -8A 



tanh 



(V^^+SA-m/cC/M-^y)) -8A 

"(^-'"I/w+tftI) 



-16 



r(r) = 
tan-type A: 

16 

^8A- (y^W+8A-mK(f(T)- jf^))' 



tanh 



-i 



VK 2 _ff 2 +8A-m K (/(T)-— ^y)J -8A 



k( H—m 



i^ + W) 



r > , 

r < , 

(115) 

r > , 

r < , 
(116) 



r(T) = < 



-1(5 



tan 



tan * 



^(ml/M+^yl-g) 



^-(v^^+SA-mKC/M-^)) 2 



+ niT 



WsA- ^Vis 2 H 2 +8A-nMs(/(r)- 7 ^y) 

tan-type B: 



fr>+7fc 



-H 



1(5 



r r = < 



^8A- (v^2l?2+8A-m K (/(r)- t^)) ' 



-16 



8A- I v / K 2 -f/ 2 +8A-m K (/(r)- -JM 



tan x 



8A-I VK 2 i/ 2 +8A-mK(/(r)-T^) 



8A-(VK 2 J f/ 2 +8A-m K (/(r)- 7 i 7 ))' 



\ 



+ 717T 



r > 



r < 



/ 

(117) 



(.(ff-ml/M+T^jl 



+ n7r 



tan 



V 



8A- [vW^+lA-rnKf/M- — Ly)"j 



H~m 



/»+7^ 



+ n7r 



r > 



r < 



/ 

(118) 

The exact r(r) solutions corresponding to the phase space trajectories in Fig. 4 are given 
by tanh-type A solution (|115|) and are plotted in Fig.6. The motions are bounded and 



21 



periodic. Comparison of three curves in Fig.6 indicates that a negative cosmological constant 
A < acts effectively as an attractive force: for the same value of H Q , the particles do not 
achieve as wide a proper separation, and the frequency of oscilliation is more rapid. As well, 
a positive A > acts as a repulsive force: the frequency of oscillation decreases and the 
particles achieve a wider proper separation. 



Fig.6 

The exact r vs t curves corresponding to the phase space trajectories in Fig. 4- 



The period T for the bounded motion is obtained from tanh-type A solution with the 
condition r = and p = ± pO: 

16 tanh -l ( ^^T8Av^^? \ q 

v / K 2 m 2 +8A V K,Hm J lm ' 



16V H 2 - 4m 2 
nHm 



lm = 0, (H9) 



M .tan" 1 ( v^CTvg S) 7m < . 



, V-K 2 m 2 -8A \ nHm 

In figure 7 we plot r{r) for fixed A = —1.5 and Hq = 16 for several different values 
of m. Though the attractive effect of a negative A is common in all cases, a special (and 
rather surprising) situation arises. As the motion becomes more relativistic (i.e. m gets 
smaller) we find that a second maximum develops in the curve (see m = 0.05 curve). The 
description of the motion is as follows. The two particles start at r = depart in opposite 
directions, reaching a maximum separation. They then go back toward one another for a 
certain period of proper time. However at some point they each reverse direction, reaching 
a second maximal separation. They then reverse direction again, finally returning to their 
starting point where the motion then repeats itself. 

As the mass becomes very small, the second maximum prevails. This peculiar behavior 
H is due to a subtle interplay between the gravitational attraction, cosmological constant 
and relativistic motion of the particles. To our knowledge it has never been previously 
observed. The changes of the peaks are clearly grasped in the phas e space trajectories in 
Fig. 8. (In Fig. 7 the first maximum of m = 0.001 curve could not be drawn due to plotting 
precision.) 



Fig.7 



A sequence of equal mass curves for A = —1.5 and Hq = 16. 
Note the presence of the second maximum for m = 0.05. 
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Fig. 8 
Change of peaks for A = —1.5 and H Q = 16 as m gets smaller. 



As the negative value of A approaches its lower bound of —k 2 H 2 /8, the form of the phase 
space trajectories transforms from an S'-shaped curve to a double peaked one and then to 
a diamond shape. Figure 9 illustrates these characteristics for the case of m = 0.5 and 
H = 100, in which a double peak structure appears for A = —70. For each trajectory in 
Fig. 9 the corresponding r(r) plot is shown in Fig. 10. 



Fig. 9 

Phase space trajectories for m = 0.5, H = 100 and different values of negative A. 



Fig. 10 

The r vs r curves corresponding the phase space trajectories in Fig. 9. 



The double-peak structure shown in Figs. 7- 10 is a consequence of having a momentum- 
dependent potential. We can gain some insight into this behaviour by computing a pertur- 
bative Hamiltonian. The structure of the determining equations (|9~3| ) and ([54]) suggests that 
we can carry out a 2-parameter expansion of the Hamiltonian in terms of k and A/k 2 . To 
the third order the result is 



H = H (p,r) + H A (p,r) 



K 



2 



= 2\lp 2 + m 2 H — (yp 2 + rn 2 — ep) 2 \r\ + -r^iyP 2 + ^ 2 — ep) 3 ?" 2 

7^ 3 , n— — 2 _ ~n4 i is _ A e P i i - A e P m2 2 

6x4 3lVP m In 2k' vV + ™ 2 16 'p 2 + m 2T 

A 2 ev 

+ T^ ■ 7 2 I 2 W2 l r l + ■ • ■ ■ ( 12 °) 

4/t d (p z + m z ) 6 i z 

It is straightforward to show that the terms in H\ have the form 

v Sl 

1 W\ S3 , (121) 



(y/p 2 + m 2 ) s 2 
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to arbitrary order in A/k 2 , where the s, are positive integers and S2 > si. One of the 
characteristics common to all such terms is that they vanish as p — > 0, since in this case the 
determining equation of the Hamiltonian becomes 

tanh (jg#Mj = — ^— • (122) 

which is A- independent. Another important characteristic is that they have a single maxi- 
mum at p 2 = Sim if So > Si in the case of A < 0. 

Consider the situation depicted in Figs. 7-10, that of two particles initially at the origin, 
each having initial momentum p , with A < 0. In this case H = 1\/p 2 + m 2 , H\ = 0, and 
the particles initially move apart as though they were free. From ( |109| ) the momentum is 
a monotonically decreasing function of time. For po/m sufficiently small, the particles will 
execute a motion which is a perturbation of that described in section |4] since the momentum 
never becomes large enough to cross the maximum in H\. However if po/m is sufficiently 
large, the terms in H\ will grow as p decreases, and the A-dependent part of the potential 
will continue to increase. The terms in Hq will decrease as p decreases, even though r is 
increasing. Eventually a maximum in r is reached, after which both r and p are decreasing. 
In the generic case r will continue to decrease toward zero. However a second extremum 
will appear if the H\ terms get too small too rapidly before p — 0, which can happen for 
A within a certain range. Since H = Hq + H\ is a constant of the motion, the only way 
to preserve the constancy of H will be for r to increase again. Essentially the particles are 
repelled due to their kinetic energy within this range of A. Of course r cannot increase too 
much, because p continues to get small - eventually r must reach a 2nd maximum, and then 
turns around again until p = 0, after which the motion continues to r = whereupon the 
particles interchange roles. 

We can see this effect in a simple non-relativistic model with 

H = p 2 /m + Km —. — -\r\ = m \ p 2 + K — It"! 1 , (123) 

1 + jr jvnr \ 1 + p J 

where in the latter equation p and r have been rescaled in units of m. The potential has 
been chosen so that S2 = 1 + s% = 2 in terms of (|121|) above. Since H = mh is a constant of 
the motion, we can write 

, (-p 2 + h)(l+p 2 ) 



m 



Kp 



(124) 

which has extrema at p — p± = | y — 6 + 6 h ± 6 a/1 — \Ah + h 2 . For p\ = h < 7 + 4v3, 
the extreme values of \r\ are not real, and so there will be no double peak structure in the 
phase-space trajectory. However for p\ = h > 7 + 4\/3 two extrema appear, and the particle 
gets a bounce. This motion is viewed qualitatively in a potential diagram shown in Fig.ll. 
As the momentum decreases from a initial value po to zero, the potential curve (which is 
linear in \r\) changes from AtoB— > ■ ■ ■ —> E — >F^E^---^B^>A^>G— ►!—»••••. 
Accordingly the particle moves from points to 10 in numerical order. 

The equations of motion for the Hamiltonian (|123|) are easily solved, yielding 

, lh - W(h e h - 2Kt )\ (l + W(he h - 2Kt )\ 

p{t) = ^W(he^) , r(t) = I ^ J \ V 1*- , (125) 

K^W{he h ~ 2Kt ) 
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where W(x) is the Lambert-W function ([SOT). The function p(t) is monotonically decreasing. 
Provided h > 7 + 4^/3 the particle gets a bounce at r = r{p_) before moving out to infinity, 
as illustrated in Figs. 12 and 13 for the case of h — (7 + 4\/3) + 20. In these figures the 
numbers on the curves denote the corresponding numbers in the potential diagram of Fig. 11. 



Fig. 11 



A schematic view of the motion with a large po in a potential diagram. 



Fig. 12 



The p(t) and r(t) curves in the non-relativistic model for h = (7 + 4-\/3) + 20. 

A bounce occurs at t = £_ = 18.08. 



Fig. 13 

The phase space trajectory in the non-relativistic model for h = (7 + 4v^3) + 20. 



The expansion (|120|) has similar features, except that there is an additional gravitational 



and cosmological attraction which prevents the separation from diverging. To order A the 
potential is a sum of two terms of the form in (|121|) , with s 2 = Si = 1 in the first term 
and s 2 = 1 + Si = 2 in the second term. The latter provides the bounce effect described 
above, but is always overwhelmed by the first term for small \r\. However the A 2 term has 
S2 = 2 + si = 3 and is a pure bounce term. Hence there exists a range of A which can provide 
a bounce. Fig. 14 shows the phase space trajectories in r > region for A = —0.3, H q = 5 and 
two different masses (m = 0.2, 0.7), in which the dotted curves represent the motions to the 
first order of A in the Hamiltonian ( |120|) and the solid curves do the motions corresponding 
to the second order of A. We see that to order A, as the motion becomes relativistic (m 
gets smaller) the trajectory simply expands and changes from S"-shape to a diamond shape. 
When the A 2 term is included, the double peak structure (a solid curve with m = 0.2) 
appears. 



Fig. 14 
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Phase space trajectories for the perturbative Hamiltonian 
for A = -0.3, #o = 5 and m = 0.2, 0.7. 



This perturbative analysis indicates that the bounce effect is a result of the negative 
cosmological constant inducing a momentum dependent potential with positive coefficients. 
For A > 0, odd powers of A/k 2 are strongly repulsive, and suppress the attractive effects of 
even powers of A/k 2 , eliminating the double peak structure. 

More generally, a positive cosmological constant acts effectively as a repulsive force. 
Figure 15 shows r(r) plots for fixed A = 1.5, Ho = 16 and several different values of m. 
The motion becomes unbounded between m = 4.72 and m = 4.73. The r(r) plots in Fig. 16 
are for fixed m = 0.5, Hq = 16 and different values of A, showing also the transition from 
bounded to unbounded motion. This transition occurs at J\ = and the critical value of A 
is given by A c = 2(H f^r 

As A — ► A c the particles rapidly separate, remaining nearly stationary for an increasingly 
large period of proper time before coming together again. At A = A c this separation time 
becomes infinite, and for A > A c , the separation diverges at finite r. 



Fig. 15 

A sequence of curves of equal mass for A = 1.5, Hq = 16. 



Fig. 16 

A sequence of curves near A c = 20.008333 for m = 7, H = 16. 



Though all the above solutions are derived from tanh-type equations (|96|) and fl97|), for a 
positive cosmological constant there exist also a countably infinite set of unbounded motions 
specified by tan-type A, B equations ( J98|) and fl99[) . Then for < A < A c , both bounded and 
unbounded motions are realized for a fixed value of H, as shown in Fig. 17. In the unbounded 
motion two particles simply approach one another at some minimal value of \r\ and then 
reverse direction toward infinity. In the trajectories the dotted curves come from tan-type 
A solution (|117|) and the dashed curves do from tan-type B solution (|118|) . As A approaches 



A c , two bulges of the solid curve (tanh-type A) and the dotted curve (tan-type A: n=0) come 
close and contact. When A exceeds A c two curves switch to the unbounded trajectories as 
shown in the solid curves in Fig. 18. The particles cross one another before receding toward 
infinity. The upper solid curve represents the motion in which p approaches the asymptotic 
values p± = ±(±y/K?H 2 + 8A + v 7 ^) as r — ► ±oo. 



26 



As noted previously, one peculiar feature of this motion is that the two particles diverges 
to infinite separation at finite proper time. The time r M for r — > oo is 



/ 



T n 



K>mJ^ m 



:log 



if (1 + J^m) - (p+ + Jp% + m 2 )(l + ^u) 



\n 



if (1 + J^) - (p + + Jp\ + m 2 )(J%^ - 1 



(126) 



The lower solid curve represents the motion in reversed direction. For A > A c only un- 
bounded motions are realized. 



Fie. 17 



Phase space trajectories of the bouned and the unbounded motions 



for A = l,?7i = 1 and if = 2.1. 



Fig. 18 

Phase space trajectories of the unbounded motions 
for A = 1.5, m = 1 and Hq = 2.1. 



We discuss in Appendix C the causal relationship between the two particles in the un- 
bounded case. 



6 THE UNEQUAL MASS CASE 



For the unequal masses the proper time ([83]) of each particle is 

16YK K imi 



dr\ 



dTo 



dt 



dt 



JKM l ^p 2 + m\ 
16YK K 2 m 2 



(127) 



JKM 2 Jp 2 + m\ 



where K = K + = K_ and Y = Y + = Y-. In this situation, choosing the time coordinate 
to be the proper time of one particle introduces an asymmetry into the description of the 
motion. Instead we seek a time variable which is symmetric with respect to 1 ^-> 2 and 
reduces to the proper time (|103|) when mi = 7772. From (127) we choose 



dr = dt 



WYKn 



KiK2mim2 



1/2 



JK \M 1 M 2 Jp 2 + m\Jp 2 + m\ 



128) 
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In terms of this variable the canonical equations are expressed as 

?. An2_|_ ™2 S 



dp l_ t K 1 K 2 M 1 M 2 ^p 2 + m 2 ^p 2 + ml 

df Ak \ m\m2 



u+l 



dzi , ^.^M 1 M 2 Jp 2 + mjJp 2 + m 2 2 \ t J K 



1/2 



+ 



V 




dr \ K\K2V(i\m2 J 1 16-Ko Mi \ l^z _)_ m ? .ft' 

(130) 



rfr / Afi M 2 ^/p 2 + m? y/p 2 + m\ 

df \ KiK 2 mim2 



1/2 



8K M x 



ip z + m{ 




W±L + #LL^-^|^|^^- fF ]|. (131) 



Note that r still describes the proper distance between the particles at any fixed instant. 



Unlike the equal mass case, the integration / dp[K\K2M\M2\i p 2 + mfyp 2 + ra\ 
can not be performed within the framework of elementary calculus. Hence we solve ( 129 ) 
numerically. 

In the case of a negative cosmological constant the rir) plots in Fig. 19 show the trajec- 
tories for various mass ratios Vfi\jvfi2 in the fixed A = — l,m 2 = 1 and H = 10. Compared 
with the equal mass case ml = ra2 = 1, as the mass ratio gets larger, the gravitational 
attraction is stronger and the proper distance between two particles as well as the period 
become shorter. When the mass ratio gets a small value than unity, the gravity becomes 
weak. However, for quite a small mass ratio a strong attractive effect of the cosmological 
constant prevails and the period changes to become shorter. At the same time the double 
peak structure (the second maximum) appears and finally the first maximum fades out. 
These characteristics are very clear in the unequal mass case. 



Fig. 19 

r{r) plots for the different values of the mass ratio m 1 /rri2 
for A = —l,m2 — 1 one? H Q = 10. 



For a positive cosmological constant the situation is simple. As shown in Fig. 20, as 
the mass ratio becomes small, the particles separate with an increasingly larger period of 
bounded motion. This is due to a repulsive effect of the cosmological constant and beyond 
the critical value the motion bocomes unbounded. 



Fig. 20 
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r(r) plots for the different values of the mass ratio m 1 /?ri 2 
for A = 0.003, m 2 = 1 and H = 10. 



7 CONCLUSIONS 

In general relativity the relationship between the motion of a set of iV bodies and the 
structure of space-time is non-linear and quite complicated, even for N = 2. Expanding 
upon the solution presented in 0, we have obtained an exact solution to the 2-body problem 
in (1 + 1) dimensions with a cosmological constant. To our knowledge this is the first non- 
perturbative relativistic curved-spacetime treatment of this problem, providing new avenues 
for investigation of one-dimensional self-gravitating systems. 
We recapitulate the main results of our paper: 

(1) We formulated the canonical formalism for a system of N bodies in a lineal theory of 
gravity with a cosmological constant A. The system is described by a conservative Hamilto- 
nian. The effect of A is incorporated into the potential. 

(2) For N = 2 the determining equation of the Hamiltonian is a transcendental equation 
derived from the matching conditions and appropriate boundary conditions at infinity. From 
these the canonical equations of motion may be derived. The metric components are also 
completely determined. 

(3) For the equal mass case we obtained explicitly the exact solutions to the canonical equa- 
tions in terms of the mutual proper time of the particles. Using the solutions we analyzed 
the motion in both r{r) plots and phase-space trajectories. 

(4) As expected, a positive cosmological constant yields a repulsive effect on the motion rela- 
tive to their mutual gravitational attraction. For < A < A c both bounded and unbounded 
motions are realized, while for A c < A only the unbounded motions are allowed. As A — > A c 
the particles separate to an infinite proper distance in infinite proper time. For A > A c this 
infinite separation occurs in finite proper time. 

(5) A negative cosmological constant has an additional attractive effect, and the motion of 
the particles is bounded. However for a certain range of the parameters, a repulsive effect 
sets in, resulting the double-peaked structures of Figs. 7- 10. This effect is due to a subtle 
interplay between the momentum-dependent A potential and the gravitational attraction. 

(6) In the unequal mass case the same basic features also occur; indeed the double peak 
behavior shows up more clearly than in the equal mass case. Although eq. ( |12ij| ) cannot be 
integrated in terms of elementary functions, it is straightforward to numerically integrate. 
An exact solution in the small mass limit of the particle 1 was also obtained. 

Several interesting features of the motion remain to be explored. The divergent separation 
of the bodies at finite proper time needs to be better understood. Another issue concerns 
the condition ( |108|) which means that for a given value of A = — |A| the motion is allowed for 



the total energy larger than J8\A\/k 2 . What is the physical meaning of this condition? It 
seems to suggest that as the attractive effect of A < exceeds a critical value the two particle 
system is no longer stable and transforms into some other system (probably making a black 
hole). To formulate the canonical formalism to treat this problem is our next subject. 
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APPENDIX A: SOLUTION OF THE METRIC TEN- 
SOR 

Under the coordinate conditions (p0|) the field equations ([12]), QT3]), ( |i~6l ) and (|TT| ) become 



1 2 8/t V ; 4 W 


~ 7= =S(X- Zl (t))- 

2yp| + m( 


P2 Sf 

2yjpb + m 2 2 


-*»(*))} 


+iVi {vr' + Pl 8{x - zx{t)) + p 2 £(z - 


z 2 (t))} + ^N^' + N[n = 


= 0, 


(132) 


KnN + N[ = , 






(133) 


9 1 (IjVo*' + J^) = 0, 






(134) 


tf + 2«iV~o7r - JVi*' = . 






(135) 


The solution to (|134j) is 













[Ml 


(+) region 


N = A e~^ 


= Acf ) 2 = 


Ml 


(0) region 






\A<f_ 


(-) region 


A being an integration constant. 








Eq. (|133|) is 









'^2 



N[ = -kA X 
The solution in each region is 

' N 1W = e{A £# - D + ) 
N m = e {A%L [Al e K °* - B 2 e~ K °*\ + 2AY A B x + D } 

where D + ,D^ and Dq are integration constants. 

The matching conditions Nx(+)(zx) = Nx(p)(zi) and Nx(-)(z 2 ) = Nx(o)(z 2 ) lead to 



{+) region 
(+) region 
[+) region 



A 
Do 



8K {D + + D_) i {KoiZl . Ko2Z2) 
J 



D+-D. 



Y+\ 1 



Yo_ _ 

K K+) Mi 



K K+) 
K K 



K x - 



M, 



Ko K. 

Yo 



Ko 



K X K 2 - -^K 1 K 2 (z 1 + z 2 ) 



(136) 



(137) 



(138) 



(139) 



(140) 



In deriving these relations the expressions ( |64"D for A +>0 , B_ and the determining equation 
(|6~6"D were used. As for the equation (|132|) , first take the 5 function at x = z\\ 



{5 function part of LHS (132) at x = z{\ 



-px 8{x - zx) < zx - N (zx) 
~pxS(x-zx)elj±-D+ 



Pi 



p\ + ml 



+ Nx(zx) 



:i41) 
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where Nq(zx), Ni(zi) and the canonical equation were inserted. Then the integration con- 
stant D + should be 

D± = l±. , (142) 



and similary 



L>_ 






(143) 



Now the metric tensor is completely determined: 
No(+)(x) 



Y± + Xz_\ K ° Kl e K+(x-zi) 



N 0H (x) 



J \K + K_J Mi 
1 fY, Y_ N 



2K J \K+ K. 
J\K + K. M 2 



(KiMi)*e 2 



1 ±k {x-zx) _|_ ^K 2 M 2 )^e^ Ko( ~ x ~ Z2) 



y + ■ Y-\KqK 2 e _ K _ {x _ Z2) 



(144) 



Y + \% (Y + , F_ ^ 1^ eK+{x _ Zi) _ i 



*«+) = ^ -F h^ + T- 



#1(0) 



#i(- 



K + \J \K + K-) Mi 

+2K (K 1 K 2 M l M 2 ) 1/2 e ^ Koizi ~ Z2) x] + A)} 



K2M2 e^ o(x " Z2) - KiMi e - Ko( ^ Zl) 



e K \j\K + + k) M 2 ' 



)l\ K K 2 c - K _ {x - Z2 ) _ 1 



With this solution and the canonical equations, the field equation ( |132| ) can be proved to 
hold in a whole x space. 

As we showed in the previous paper, to satisfy (|135|) the dilaton field \l/ needs an extra 



function /(£), which has no effect on the dynamics of particles. After lengthy calculation 
Eq.(PS|) leads to 



Pi 



2K n K- 



0^2 



+ 4 ( AoA , + KoK2 _ ftftft _ ftftft) + 4cYoKd (£ + ft 



Thus /(£) is uniquely determined. 



P2 



\jvl + m 2 2 



(145) 
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APPENDIX B: A TEST PARTICLE APPROXIMA- 
TION 

For a single static source M the solution to the field equations (^) -fllTD under the coordinate 
conditions fl20T) is 



No 



7[ 



,^N 



N x = e\l 



8A a; 



k 2 M 2 Id 



,fN _ ! 



(146) 



*, + 8A 



K' 



M 2 



kM . , /teM / 8A 



{}From ([Tj^) and (|TU|) the canonical equations for a test particle (mass //) under the gravity 
of a static source are 



P 



-Vp +/ 1 



6*2 



+ p 



6*2 



p 



Vp 2 + /i 2 

The Hamiltonian leading to these equations is 



= N (z)-N 1 (z) 



H = M+Jp 2 + fi 2 N (z)-pN 1 (z) 



= M + Jp 2 + /x 2 e— |z| - em/1 + -r^ -- e~ |z| - 1 
v V fir M" 1 |z| L 

This Hamiltonian is also derived from the determining eq. ( |5B] ) by setting 

z 
Z! = z, m 1 =n, pi=p, Pi=p=p-r-r, 

\z\ 

z 2 = 0, m 2 = M, p 2 = 0, 
and retaining only the linear terms of y/p 2 + \x 2 and p. 

In terms of the proper time of the test particle 



dr 2 = dt 2 {n (z) 2 - (Nx(z) + zf) = dt 2 N ( 



P 



°^2 



P 1 + jJL 2 



the canonical equations ( |147| ) and ( |148| ) are expressed as 
dp kM 



dr 

dz 
dr~ 



Afi 



V + H \ vP + A* sgn(z) - epy 1 + 



8A 

K 2 M 2 



P 



„ i l+ ^M 2 ^ z) ^- e * w )— ^- 



(147) 
(148) 



(149) 



(150) 



:i5r 



(152) 
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Eq. (|151|) can be integrated and in z > region the solution p(r) is 

p(r) 



?K*5 



with 



h(r) 



(l+VW) l-peVft^' 



v 7 ^ 



Til+p e e V^ (T - T 0) 






i- (1+ ^ I7)( V^!l^tan 



(t-to) 



y^I+g 



A>0 



A = 



A<0 



lf^ + ^^ tan 



eV-f^-^o) 



where 



8A 



P 



(1 + v / 7m)(v // Po + ^ 2 ~ epo) ~ /V7m - 1 
(1 + v^Tm) (\Ipo + /^ 2 - e Po) + /V7m - 1 



with po being the initial momentum at r = Tq. 
In z < region the solution is 



P(r) = -f[h(r) 



h(r) 



with 



h{r) 



{ 1 +\fTM){ l—pe 



#(--nj) 



V7A/-1S 1+P e 



4< t - t o) 






1 (i+v^uXy/pg+M^.Po) tan 



e-v/^^-l)) 



■ v /p§+M 2 +€ P0 0^77 



A>0 



A = 



A<0 



i+V^M 



tan 



|(r-ro) 



where 



1 + \/Jm)(Jpq + A* 2 + epo) - HV7M ~ 1 



(153) 



(154) 



(155) 



(156) 



(157) 



:i58i 



(1 + Vim){\JpI + /i 2 + epo) + nVlM - 1 
When p = Po and z = at r = r , the total energy is H = H = M + J pi + /i 2 . The 
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solution for z(t) is obtained from (|149|) and p(r) as 



z T 



Vpq+m 2 



-(fe-i)VT^" 



k m log (h+i) _ (M)vt5? 



log 



Vpg+M 2 



(*-£)VW 



«M -6 (^i)-(fe-i)v^F 



For the test particle solution the critical value of A is A c 



K fl' 



z>0 



z <0 . 



■M 2 



(159) 



H ■ 



Fig.21 and 22 show typical trajectories of the test particle fi = 0.1 for M — 10 and 
A = —10, 0, 0.4 and 2. The characteristics of these plots are common to those of the unequal 
mass case. 



Fig.21 



Phase space trajectories of a test particle for different values of A. 



Fig. 22 



The r plots corresponding to the trajectories in Fig.21. 



APPENDIX C: CAUSAL RELATIONSHIPS BETWEEN 
PARTICLES IN UNBOUNDED MOTION 

We can explicitly verify that the particles lose causal contact with one another for r > t^. 
Consider the unbounded motion of tan-type A(n = 0) with m — 1,H — 2.1 and A = 1.5. 
The path x(t) of light emitted from particle 2 at time T is governed by dr = 0, which reads 



(I) + w "i - w - n >> - ° 

and so the equation of the light signal directed to particle 1 is 

dx 

— = N {x{t),z 1 {t),z 2 {t),p{t))-N 1 {x{t),z 1 {t),z 2 {t),p{t)) 



(160) 



(161) 
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The light emitted in the opposite direction is described by 

dx 

— = -N (x(t),z l (t),z 2 (t),p(t)) - N 1 (x(t),z l (t),z 2 (t),p(t)) . (162) 

Numerically solving ( |161|) and (|162|) yields the solutions shown in Figs. 23 and 24, where the 
trajectories of light signals emitted from particle 2 at various times T are plotted. For small 
T(T < 3), the particles are in causal contact (a dotted curve in (+) direction in Fig. 23), 
but for T p« 3 the signal just barely catches up with particle 1, which is almost in light-like 
motion (a dashed curve in (+) direction in Fig.23). For T = 4 the world line x(t) in the 
(+) direction is parallel to Z\(t) at large t and in the (-) direction it goes nearly on the same 
trajectory with the particle 2. For large T (T > 4.82) the particles are out of causal contact 
with each other (Fig.24): a light ray sent from particle 2 toward particle 1 receives a strong 
repulsive effect and ultimately reverses direction, following behind particle 2. In Fig.24 the 
trajectories of the light signal emitted to (-) direction can not be discriminated from those 
of the particle 2. 



Fig.23 

The trajectories of light signals emitted at T = 2, 3 and 4. 



Fig.24 
The trajectories of light signals emitted at T = 4.82, 5 and 6. 



A flat-space model of these effects can be constructed as follows. Consider the following 
expression for the 2- velocity 

«" = (/Kl//V)-l) (163) 

where /(err) is some function and 

dr 2 = dt 2 - dx 2 (164) 



is the flat metric. We have j- = /, ¥- = V ' f 2 — 1 and 



so 



Tt = — /— ' (165) 

The general expression for the acceleration of a particle with 2- velocity ( |163j ) is 

* f{l, -=£==) (166) 



a 



A' 



dr v// 2 (ar) - 1 
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where /' = df(r)/dr. We have u ■ u = 1 and a ■ u = and 

{Off , 1R7 x 

ffl - ffl= /V)-i (16T) 

for the magnitude of the acceleration. In general we have the following possibilites: 

1) The function / — > f where /o is finite at r — > oo. In this case the particle never 
becomes lightlike. 

2) The function / — > oo as r — > oo. In this case the particle becomes lightlike, but it takes 
an infinite proper time (and coordinate time) for this to happen. The standard example is 
/ = cosh (err), the constant acceleration example. 

3) The function / — > oo as r — > r , where r is finite. In this case the particle becomes 
lightlike in a finite amount of proper time, but an infinite amount of coordinate time. An 
example would be / = sec(crr). The acceleration is not constant, but increases as a function 
of proper time, diverging at r = To- This last situation is realized by our exact solutions 

iT7HTT8|) with A > A c . 
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